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ABSTRACT 

We present the current photometric dataset for the Sloan Lens ACS (SLACS) Survey, including 
HST photometry from ACS, WFPC2, and NICMOS. These data have enabled the confirmation of an 
additional 15 grade 'A' (certain) lens systems, bringing the number of SLACS grade 'A' lenses to 85; 
including 13 grade 'B' (likely) systems, SLACS has identified nearly 100 lenses and lens candidates. 
Approximately 80% of the grade 'A' systems have elliptical morphologies while ~10% show spiral 
structure; the remaining lenses have lenticular morphologies. Spectroscopic redshifts for the lens and 
source are available for every system, making SLACS the largest homogeneous dataset of galaxy-scale 
lenses to date. We have created lens models using singular isothermal ellipsoid mass distributions for 
the 11 new systems that are dominated by a single mass component and where the multiple images 
are detected with sufficient signal-to-noise; these models give a high precision measurement of the 
mass within the Einstein radius of each lens. We have developed a novel Bayesian stellar population 
analysis code to determine robust stellar masses with accurate error estimates. We apply this code 
to deep, high-resolution HST imaging and determine stellar masses with typical statistical errors of 
0.1 dex; we find that these stellar masses are unbiased compared to estimates obtained using SDSS 
photometry, provided that informative priors are used. The stellar masses range from 10 10 5 to 10 118 
Mq and the typical stellar mass fraction within the Einstein radius is 0.4, assuming a Chabrier IMF. 
The ensemble properties of the SLACS lens galaxies, e.g. stellar masses and projected ellipticities, 
appear to be indistinguishable from other SDSS galaxies with similar stellar velocity dispersions. This 
further supports that SLACS lenses are representative of the overall population of massive early-type 
galaxies with M* > 10 11 Mq, and are therefore an ideal dataset to investigate the kpc-scale distribution 
of luminous and dark matter in galaxies out to z ~ 0.5. 

Subject headings: galaxies: elliptical and lenticular, cD - gravitational lensing - surveys 



1. INTRODUCTION 

Early-type galaxies have been found to be remarkably 
homogeneous, as evidenced by the ti ght scaling relation- 
ships of the Fundamental Plane (jDjorgovski fc Davisl 
Il987t iDressler et al.l I1987D and the small scatter 
of the ir nearly-isotherm al cen tr al density p rofiles 
(e.g.. iRusin fe Kochanekl 120051; iRusin et al.l 120031 ; 
iKoopmans et all I2006L 12009) . However, there are many 
details of this homogeneity that remain unresolved, 
including an explanation for the tilt of the Fundamental 
Plane with respect to the virial plane and f or the scatter, 
or fini te thickness, orthogonal to the plane. iBolton et al.l 
( 2008) have used lensing and stellar velocity dispersions 
to show that the dynamical mass scales with the lensing 
(i.e., total) mass of early- type galaxies, suggesting that 
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non- homology is unlikely to cause the observed tilt. 
Joint dynamical and stellar population modeling has 
suggested that more massive earl y- type galaxies have 
highe r dark matter fractions (e.g., iPadmanabhan et al 



2004) that could explain most of the tilt (jTortora et al 
2009). These models depend on assumptions about 



the initial mass function and do not readily explain 
how the luminous bulge and a dark matter halo might 
'conspire' to produce an isothermal total mass density 
profile; a well-defined set of galaxies with accurate 
photometric data is required in order to perform the 
detailed modeling necessary to accurately separate the 
dark and luminous mass components. 

A particularly powerful tool is the combination 
of precise strong gravitational lensing mass measure- 
ments, stellar velocity dispersions, and accurate stellar 
mass estimates. The advent of the Sloan Lens ACS 
Survey (SLACS: IBolton et all l200l ITreu et all l2(M 



I Koopmans et al.ll2006l; IGavazzi et alJl2007t IBolton et al 
2008[ iGavazzi et al.ll2008t IBolton et alj|2008t ITreu et al 
120091 papers I- VIII, respectively) has made this type 
of modeling feasible for the first time for a large and 
uniformly selected dataset. One of the key features of 
SLACS is that each lens contains all of the information 
necessary to perform a robust decomposition of the mat- 
ter content in the central regions of galaxies. For ex- 
ample, the SDSS spectroscopy provides a precise mea- 
surement of the stellar velocity dispersion and the lens 
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and source redshifts, while high-resolution Hubble Space 
Telescope (HST) imaging yields a detailed view of the 
lensed background source and the surface brightness pro- 
file of the lensing galaxy. 

In this paper we present the definitive SLACS sam- 
ple and an essential component necessary for disentan- 
gling the dark and luminous matter in early-type galax- 
ies: multi-band optical and near-infrared imaging from 
HST. The photometric colors provided by HST allow us 
to accurately constrain the stellar mass of the lensing 
galaxy, thereby determining the normalization of the lu- 
minous component of the mass distribution. Although 
other authors have investigated the stellar populations 
of SLACS lenses using SDSS photometry (|Grillo et al.1 
120081 120091 hereafter G09), we find that the HST imag- 
ing is necessary to adequately perform model photometry 
and yield unbiased stellar mass estimates, to the level of 
accuracy that is comparable with that of the other ob- 
servables. 

The paper is organized as follows. We begin by de- 
scribing our new HST observations in Section 2. We 
then provide a detailed summary of the SLACS dataset 
in Section 3, including the introduction of 15 newly con- 
firmed strong lenses. In Section 4 we discuss our model 
photometry, and our stellar mass estimates are detailed 
in Section 5. A discussion of the ensemble properties of 
SLACS lenses, including comparison with 'twin' galaxies 
from the SDSS, is given in Section 6 before concluding 
with a summary of the final SLACS d ataset in S ection 7. 
All magnitudes are on the AB scale (jOkejll974l L a stan- 
dard concordance cosmology with Q m = 0.3 &A = 0.7 
h = 0.7 is assumed, and base-10 logarithms are denoted 
by 'log' while natural logs are written as 'In.' 

2. MULTI-COLOR IMAGING 

A campaign was initiated to obtain HST multi-band 
imaging in V, I, and H (F555W, F814W, and F160W) 
for all of the confirmed SLACS lenses and a small num- 
ber of high-probability lens candidates from Paper V us- 
ing the ACS WFC and NICMOS NIC2 cameras (Pro- 
grams 10494, 10798, and 11202; PI Koopmans). The 
programs were switched to WFPC2 after the failure of 
ACS, and F555W was replaced with F606W to exploit 
the improved throughput in F606W compared to F555W 
on WFPC2; we refer to these two filters interchangeably 
as the V-band throughout the rest of this paper. Further- 
more, the failure of NICMOS in September 2008 has left 
29 lenses unobserved in the H band. All of the ACS data 
were reduced using the methods described in Paper V. 
A total of 77 systems were observed with WFPC2 and 
59 systems were observed with NICMOS; this includes 
12 new lens candidates not included in Paper V (Table 
[T]) . A summary of the imaging for our targets is given in 
Table H 

2.1. WFPC2 Data Reduction 

The SLACS WFPC2 images were reduced using a 
custom-built pipeline that largely follows the procedures 
used in the PyRAF task MULTIDRIZZLE. One significant 
difference is the initial step; we run a Python implemen- 
tation of the LACOSMIC program on each of the raw sci- 
ence frames to identify and mask cosmic rays. The in- 
dividual exposures are then background subtracted and 
the CROSSDRiz task is used to find the sub-pixel shifts 



between exposures. We only use the shifts determined 
from the WFC3 chip because we know this chip contains 
a bright source (the SLACS target) to ensure robust cross 
correlation. This avoids problems with spurious offsets 
from the other chips (due to no objects in the chip frame, 
for example) skewing the average offset; this procedure 
yields image registration to better than 0.1 pixel. The 4 
exposures for each chip are drizzled to an output image 
using a 0'.'05 pixel scale and a pixfrac of 0.6. The drizzled 
images of the four chips are then placed into a common 
image using SWARP. All of the lens-galaxy photometry 
presented in this paper is performed on the pre-SWARP 
drizzled image of the WFC3 chip. Absolute astrometry 
is achieved by setting the location of the lens galaxy to 
the SDSS coordinates. 

2.2. NICMOS Data Reduction 

The SLACS lens candidates tend to be quite large 
compared to the NICMOS field of view, and as a re- 
sult special care must be taken when reducing the NIC- 
MOS imaging. We have used the most recent version 
of the CALNICA pre-processing pipeline (version 4.4.1, 
implemented in the OPUS pipeline and applied by re- 
downloading the _cal.fits images from the HST archive) 
to perform the flatfielding, de-biasing, linearity correc- 
tions, time-stream cosmic ray rejection, bad-pixel mask- 
ing, and error model creation. 

One particular problem presented by the large size 
of the SLACS galaxies is that the unknown floating 
bias of each amplifier (the 'pedestal') cannot be accu- 
rately quantified using standard techniques if the lens 
dominates the amplifier quadrant. We use an itera- 
tive approach to determine an approximation of the true 
pedestal by first using the STSDAS task pedsub to re- 
move an initial estimate of the pedestal. We next fit a 
model of the amplifier glow in the pedestal-subtracted 
images to remove residuals which were not adequately 
subtracted by CALNICA. The scaled amplifier glow 
model is then subtracted from the original data and PED- 
SUB is performed again on the improved amplifier glow- 
corrected images. At this point most of the instrumental 
signature has been removed from the data and we ap- 
ply the count-rate non-linearity correction implemented 
in the Python task RNLINCOR. A robust clipping/fitting 
algorithm is then employed to determine the sky level in 
each exposure; this sky is subtracted and the STSDAS 
CROSSDRiz routine is used to determine offsets between 
exposures, yielding registration to better than 0.1 pixel. 
The exposures are then drizzled into a median image 
using the STSDAS drizzle task. As in the standard 
reduction pipelines, the median image is blotted back 
to the original CCD frame for residual cosmic ray/bad 
pixel detection. However, we also use the blotted median 
image to subtract a model of the astrophysical sources, 
leaving residuals which are dominated by the quadrant- 
dependent errors in the pedestal removal. The residual 
bias in each quadrant is found by an aggressively clipped 
mean and removed. The median image used as a model 
for the astrophysical sources is, in practice, biased by 
the incorrect initial pedestal removal and further itera- 
tions should lead to progressively better de-biased data. 
However, we have found that one iteration is adequate to 
remove most of the residual bias. This procedure yields 
four exposures for each SLACS target which have been 
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Table 1 

New SLACS Targets 



Name Right Ascension Declination zj g Z\, g crgrjgg (km s 1 ) Grade 



SDSSJ1038- 


-0024 


10h38m58.45s 


-00° 


24m38.6s 


0.219 


1.306 


198±11 


B 


SDSSJ1306- 


h0600 


13h06ml3.65s 


06° 


00m22.1s 


0.173 


0.472 


237±16 


A 


SDSSJ1313- 


h4615 


13hl3m02.93s 


1(5° 


15ml3.6s 


0.185 


0.514 


262±18 


A 


SDSSJ1318- 


-0313 


13hl8m39.33s 


-03° 


13m34.2s 


0.240 


1.300 


213±17 


A 


SDSSJ1319- 


hi 504 


13hl9m32.50s 


15° 


04m26.0s 


0.154 


0.606 


235±13 


A 


SDSSJ1325- 


h4145 


13h25m55.31s 


41° 


45m02.2s 


0.306 


0.785 


259±15 


B 


SDSSJ1436+3640 


14h36m05.63s 


36° 


40m47.1s 


0.185 


0.758 


280±18 


A 


SDSSJ1527- 


h0314 


15h27m47.69s 


03° 


14m32.1s 


0.222 


0.538 


231±16 


X 


SDSSJ1644- 


S-2625 


16h44m43.09s 


26° 


25m25.4s 


0.137 


0.610 


229±11 


A 


SDSSJ1709- 


h2324 


17h09m38.96s 


23° 


24m08.4s 


0.347 


0.719 


286±30 


A 


SDSSJ1719- 


h2939 


17hl9m34.16s 


29° 


39m26.4s 


0.181 


0.578 


286±15 


A 


SDSSJ2343- 


-0030 


23h43m58.87s 


-00° 


30m22.4s 


0.181 


0.463 


268±10 


A 



Table 2 

SLACS Imaging 



Photometry 


# A Lenses 


# B Lenses 


# X Targets 


F435W 


33 


6 


18 


F555W 


16 








F606W 


61 


8 


5 


F814W 


82 


11 


11 


F160W 


51 


3 


1 


One band 


6 


1 


24 


Two bands 


22 


3 


19 


Three bands 


32 


6 


2 


Four bands 


25 









dark-, amplifier glow-, pedestal-, and sky background- 
subtracted, and these pointings are then drizzled on to 
a common output image with a scale of 0'.'05 using a 
pixfrac of 0.8. 

We also take care in assuring that we have a robust 
estimate of the variance in this final drizzled image. 
We convert the ERR frames created by the CALINCA 
pipeline into variance images by squaring them; these 
variance images contain an estimate of the noise prop- 
erties from the dark current and amplifier read noise. 
The single-exposure variance images are then added to 
the (non-background-subtracted, pre-drizzle) science im- 
ages, which have been appropriately scaled by the ex- 
posure time per pixel (encoded in the TIME extension 
created by CALNICA) and the detector gain g, to cre- 
ate pixel-to-pixel variance maps for each exposure. These 
variance maps are individually drizzled to the output sci- 
ence coordinate frame, yielding single-exposure drizzled 
variance maps (VAR) and the drizzle 'weight' (WHT) im- 
ages. The WHT images are then used to determine the 
final variance image associated with the drizzled science 
image by taking a weighted sum of the VAR images. We 
expect this procedure yields the optimal representation 
of the variance in each pixel of the output (stacked and 
drizzled) science image. 

3. THE DEFINITIVE SLACS SAMPLE IN COLOR 

The SLACS survey has observed 12 new targets (Table 
[T]) since the conclusion of the ACS program (Paper V). 
The WFPC2 F606W images of these targets are shown 
in Figure [U including the residuals from B-spline galaxy 
model subtractions (see Section 3.1). In addition to these 
12 candidates, we revisit 14 systems from Paper V that 
had grades poorer than 'A' for which we have obtained 
new multi-band imaging. Grade 'A' refers to confirmed 



lenses, while grade 'B' denotes probable lenses and grade 
'X' is used for unlikely lenses; see Paper V for a more 
detailed description of the grading scheme. 

The new NICMOS and WFCP2 data for these sys- 
tems were treated in a similar manner to the ACS data 
of Paper V in order to establish the lensing nature of 
each system. The putative lensing galaxy was subtracted 
from each image using B-spline models, and the residu- 
als were visually inspected for evidence of lensed images 
and counter-images. A consensus grade was compiled 
by four authors (MWA, ASB, LAM, and TT) using the 
same grading scheme as Paper V. Four of the systems 
were initially given a 'B' grade by the consensus anal- 
ysis but were subsequently upgraded to 'A' when suc- 
cessful lens models were fit to the images (see below). 
Furthermore, SDSSJ1313+0506 has been confirmed as a 
lens based upon the identification of a counter-image to 
the background [O II] emission from a DEIMOS spec- 
trum (A. Dutton, private communication). In total, 9 of 
the 12 previously unobserved candidates have been con- 
firmed as 'A' lenses, and 6 systems reported on in Paper 
V have been upgraded to 'A' status. This brings the 
number of grade 'A' lens systems identified by SLACS 
to 85 (Table [3]), with an additional 13 'B' candidates. A 
color montage of the grade 'A' lenses can be found in 
Figured 
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Figure 1. The WFPC2 F606W image of the SLACS candidate 
SDSSJ2343-0030. The left panel shows the early-type morphology 
of the lens and evidence for the lensed nature of the background 
source. The right panel shows the system with B-spline model 
of the lensing galaxy subtracted, revealing a clear four-image lens 
system with a partial Einstein ring. See the electronic version of 
the journal for images of all 12 of the new targets. 



3.1. Lensing Analysis 



4 



Auger et al. 



□ 


• 


1 






1 ■flllUU 

□ 


r J I'VniHW 





ratlin 4iqHi4ii 





# • m 

m _ ■ 




Figure 2. Color cutouts of 84 of the SLACS lenses, ordered by redshift. The lens SDSSJ1618+4353 has two primary lensing galaxies and 
has not been included. See the Journal version for full-resolution figure. 

We produce singular isothermal ellipsoid (SIE) lens 



models (e.g., IKormann etHl Il99l for 11 of the 15 
new grade A lenses. One candidate from Paper V, 
SDSSJ1029+6115, is shown clearly to be a lens using IFU 
data but our HST data do not have sufficient sensitivity 
to adequately model the lensing, nor do we model the 
spectroscopically confirmed lens SDSSJ1313+0506. Ad- 
ditionally two new candidates from our WFPC2 imag- 
ing, SDSSJ1319+1504 and SDSSJ1436+3640, will re- 
quire special subtractions due to the presence of nearby 
galaxies in the field. 

The SIE modeling is performed in the same manner 
as Paper V, with the only difference being the switch 
from ACS data to WFPC2 data. To briefly summarize, 
the mass of the lensing galaxy for each system is mod- 
eled as an SIE profile positioned on the center of the 
light distribution. A lens galaxy-subtracted image of the 
background source is created by modeling the light of the 
lens galaxy using a radial B-spline model and subtract- 
ing this from our HST images. The background source is 
modeled as single or multiple Gaussian or Sersic ellipsoid 
components which are lensed by ray-tracing through the 
SIE mass model of the lens; this lensed source is fit to 
the observed images in the lens galaxy-subtracted frame. 
The model parameters (the Einstein radius and elliptic- 
ity parameters for the lens) are first adjusted by hand 
before optimization by MPFIT (jMarkwar dt 2009). The 

9 These data will all be available online at http:/ /www. slacs.org 



resulting models are shown in Figure [3] and the param- 
eters of the lens models are listed in Table |U We note 
that these are not perfect lens models (not all of the 
lens are isothermal, nor are the sources always ideally 
described by Sersic profiles). However, the SIE models 
do provide accurate determinations of the Einstein radii 
and the masses within the Einstein radii; these are the 
parameters of interest for this paper. 

3.2. The SLACS Catalogs 

The SLACS data products can be split into two cat- 
egories 9 : 'basic' data products which include photomet- 
ric and spectroscopic measurements, and inferred data 
products, including parameters inferred from lens mod- 
eling and stellar population analysis. In Table [3] we 
present our 'basic' data for the grade 'A' lenses, includ- 
ing the location on the sky, the lens and source redshifts, 
the lensing galaxy velocity dispersion, the morphological 
type, and the effective radii (measured at the intermedi- 
ate axis) and magnitudes from our photometric analysis 
(Section 4). Table Q] presents the inferred quantities, in- 
cluding Einstein radii and masses, stellar mass estimates 
for Chabrier and Salpeter initial mass functions (IMFs) 
denoted with the subscripts 'Chab' and 'Salp' respec- 
tively, luminosities, and /* ) Ein ) the stellar fraction within 
the Einstein radius for each IMF. 
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Figure 3. The SIE lens model for the SLACS lens SDSSJ2343-0030. The five panels show: the HST image of the lens; the lensing galaxy- 
subtracted residuals; the model of the lensed source and the lensing critical curves (in white); the source- and lens-subtracted residuals; 
and the reconstructed background source and caustics (in white). See the electronic version of the journal for all 11 lens models. 

Table 3 

Photometric and Spectroscopic Properties of SLACS Lenses 



Name 


RA 




Dec 




z 




%src 


(km s" 1 ) 


Morph 


m B 


r e .B 
(") 




?e,V 
(") 


mj 


r 

( 


•,I 




SDSSJ0008-0004 


00h08m02 


96s 


-00°04m08 


2s 





440 


1.192 


193±36 


E 








19.60+ 


1.82 


18.52 


1 


71 






SDSSJ0029-0055 


00h29m07 


78s 


-00°55m50 


5s 





227 


0.931 


229±18 


E 






17.74+ 


2.58 


17.05 


2 


16 






SDSSJ0037-0942 


00h37m53 


21s 


-09°42m20 


Is 





195 


0.632 


279±10 


E 


18.57 


2.87 


16.90t 


2.68 


16.20 


1 


80 


15 


19 


SDSSJ0044+0113 


00h44m02 


90s 


01°13ml2 


6s 





120 


0.197 


266±13 


E 






16.32t 


3.25 


15.69 


1 


92 


15 


03 


SUSS J 0109+1500 


01h09m33 


73s 


i ro r\r\ o o 

15 00m32 


5 s 





294 


0.525 


251±19 


E 






" ' 




17.63 


1 


38 






SDSSJ0157— 0056 


01h57m58 


94s 


-00°56m26 


Is 





513 


0.924 


295±47 


E 






19.87 T 


1.21 


18.59 


1 


84 


17 


52 


SDSSJ0216— 0813 


02hl6m52 


54s 


-08°13m45 


3s 





332 


0.523 


333±23 


E 


19.91 


2.49 


18.36* 


2.97 


16.86 


2 


40 


15 


99 


SDSSJ0252+0039 


02h52m45 


21s 


00°39m58 


4s 





280 


0.982 


164±12 


E 






18.77+ 


1.36 


17.87 


1 


39 






SDSSJ0330— 0020 


03n30ml2 


14s 


-00°20m51 


9s 





351 


1.071 


212±21 


E 






18.84+ 


1.52 


17.99 





91 


17 
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Table 3 

Photometric and Spectroscopic Properties of SLACS Lenses 



SDSSJ1319+1504 
SDSSJ1330-0148 
SDSSJ1402+6321 
SDSSJ1403+0006 
SDSSJ1416+5136 
SDSSJ1420+6019 
SDSSJ1430+4105 
SDSSJ1432+6317 
SDSSJ1436+3640 
SDSSJ1436-0000 
SDSSJ1443+0304 
SDSSJ1451-0239 
SDSSJ1525+3327 
SDSSJ1531-0105 
SDSSJ1538+5817 
SDSSJ1614+4522 
SDSSJ1618+4353 
SDSSJ1621+3931 
SDSSJ1627-0053 
SDSSJ1630+4520 
SDSSJ1636+4707 
SDSSJ1644+2625 
SDSSJ1709+2324 
SDSSJ1718+6424 
SDSSJ1719+2939 
SDSSJ2141-0001 
SDSSJ2238-0754 
SDSSJ2300+0022 
SDSSJ2302-0840 
SDSSJ2303+1422 
SDSSJ2321-0939 
SDSSJ2341+0000 
SDSSJ2343-0030 
SDSSJ2347-0005 



13hl9m32.50s 
13h30m45.53s 
14h02m28.21s 
14h03m29.49s 
14hl6m22.34s 
14h20ml5.85s 
14h30m04.10s 
14h32ml3.34s 
14h36m05.63s 
14h36m27.54s 
14h43ml9.62s 
14h51m28.19s 
15h25m06.70s 
15h31m50.07s 
15h38ml2.92s 
16hl4m37.74s 
16hl8m43.10s 
16h21m32.99s 
16h27m46.45s 
16h30m28.16s 
16h36m02.62s 
16h44m43.09s 
17h09m38.96s 
17hl8m37.39s 
17hl9m34.16s 
21h41m54.68s 
22h38m40.20s 
23h00m53.15s 
23h02m20.18s 
23h03m21.72s 
23h21m20.93s 
23h41mll.57s 
23h43m58.87s 
23h47m28.08s 



lo U4mzb.Us 


0.154 


0.606 


235±14 


Ul 4em41.os 


0.081 


0.711 


185± 9 


n o o o 1 o o rr„ 

bo 21moo.5s 


0.205 


0.481 


267±17 


UU Uom41.4s 


0.189 


0.473 


213il7 


cr i o c on ,1 _ 

51 obmoU.4s 


0.299 


0.811 


240±25 


bU lyml4.os 


0.063 


0.535 


205i 4 


a 1 o n rr rr 1 „ 

41 UomoY.ls 


0.285 


0.575 


612±62 


n o o -i t no o 

63 17m03.8s 


0.123 


0.664 


199± 8 


36 40m47.1s 


0.185 


0.758 


281±19 


c\r\Q c\f\ nn n 

00 00m29.2s 


0.285 


0.805 


224±17 


03°04m08.2s 


0.134 


0.419 


209±11 


noo on _ on a . , 


0.125 


0.520 


223il4 


oooot /i ^ /i 

33 27m47. 4s 


0.358 


0.717 


nn a 1 n n 

264±26 


01°05m45.7s 


0.160 


0.744 


279±12 


58 17m09.8s 


0.143 


0.531 


189±12 


45°22m53.3s 


0.178 


0.811 


182±13 


43°53m27.4s 


0.199 


0.666 


292±29 


39°31m44.6s 


0.245 


0.602 


236±20 


00°53m57.6s 


0.208 


0.524 


290±14 


45°20m36.3s 


0.248 


0.793 


276±16 


47°07m29.6s 


0.228 


0.675 


231±15 


26°25m25.4s 


0.137 


0.610 


229±12 


23°24m08.4s 


0.347 


0.719 


286±30 


64°24m52.2s 


0.090 


0.737 


273±16 


29°39m26.4s 


0.181 


0.578 


286±15 


00°01ml2.3s 


0.138 


0.713 


181±14 


07°54m56.0s 


0.137 


0.713 


198±11 


00°22m38.0s 


0.228 


0.463 


279±17 


08°40m49.5s 


0.090 


0.222 


237± 8 


14°22ml7.9s 


0.155 


0.517 


255±16 


09°39ml0.3s 


0.082 


0.532 


249± 8 


00°00ml8.7s 


0.186 


0.807 


207±13 


00°30m22.4s 


0.181 


0.463 


269±11 


00°05m21.3s 


0.417 


0.714 


404±59 



E/S0 
SO 
E 
E 
E 
SO 
E 
S 

E 
SO 
E 
E 
E 
E 
E 

E 
E 
E 
E 
E 
E 

E/S0 
S 
E 
E 
E 
E 
E 
E 

E/S0 
E 



18.91 1.00 
19.07 2.56 



16.; 



2.35 



19.21 1.13 
17.88 3.49 



19.27 2.53 

19.27 2.85 

19.42 2.02 

19.76 1.81 

19.62 1.76 



17.39 4.1 



18.21 3.19 

19.80 1.86 

17.96 4.28 

16.52 5.54 



20.73 1.81 



17.89+ 
17.56+ 
17.51* 
17.64t 
18.58+ 
15.94* 
17.80* 
15.78* 
16.95* 
17.94* 
17.62* 
16.92* 
18.10* 
16.43* 
17.42* 
17.50* 

17.55* 
17.87* 
18.08* 
17.79* 
16.85* 
18.51* 
16.15* 
17.52* 

17.18* 
18.19* 
15.95* 
16.77* 
15.27* 
17.02* 
17.17* 
18.82* 



1.44 
1.36 
2.62 
1.87 
1.38 
2.12 
2.71 
6.21 
4.64 
2.94 
1.38 
2.64 
3.40 
3.32 
1.59 
2.93 

2.77 
2.05 
2.05 
1.84 
2.32 
2.06 
4.51 
1.83 

2.41 
1.93 
3.18 
3.54 
4.79 
4.52 
2.74 
1.80 



17.23 


0.77 


16.32 


16.92 


0.96 


16.18 


16.30 


2.29 


15.55 


17.02 


1.14 


16.11 


17.55 


0.98 


16.70 


15.05 


2.25 


14.37 


16.85 


2.55 




15.13 


3.04 


14.53 


17.16 


1.63 


16.38 


17.00 


0.70 


16.24 


15.92 


1.54 


15.27 


17.07 


2.42 


16.09 


15.79 


1.97 


15.14 


16.63 


1.00 


15.99 


16.80 


2.58 




17.08 


1.80 




16.79 


1.51 


16.08 


16.71 


1.98 




16.77 


1.65 


15.80 


16.98 


1.68 




16.22 


1.55 


15.46 


15.39 


1.70 


15.35 


16.87 


1.46 


15.99 


16.72 


1.81 




16.12 


1.82 


15.40 


16.96 


1.52 


16.15 


15.45 


1.70 


14.80 


15.72 


2.94 


15.05 


14.61 


4.11 




16.30 


2.36 


15.40 


17.83 


1.14 


16.85 



T From the F606W filter on WFPC2. 

* From the F555W filter on ACS. 

* The lens galaxy has at least one close companion galaxy. 

4. MODEL PHOTOMETRY 

The program GALFIT (jPeng et all I2002D is used to 
perform de Vaucouleurs model fitting to the WFPC2 and 
NICMOS images of the lenses; we adopt the model pho- 
tometry using the method from Paper V for the ACS 
imaging. Subimages of 30" on a side are used to fit the 
WFPC2 data, and we use the region of the NICMOS 
mosaic that was covered in all four pointings (generally 
an area approximately 15" x 15"). Masks of neighboring 
galaxies and the flux from any lensed images are manu- 
ally created by iteratively modeling and subtracting the 
lens galaxy. Bright galaxies very close to the primary 
lensing galaxy are simultaneously modeled with GAL- 
FIT to accurately separate the two objects. We test the 
robustness of the de Vaucouleurs model magnitudes ob- 
tained with GALFIT by comparing the results for lenses 
observed multiple times and by comparing with magni- 
tudes from a Sersic fit with the Sersic index left as a free 
parameter; we find that the de Vaucouleurs magnitudes 
are consistent with a scatter of ~ 0.03 magnitudes. We 
also record the effective radii, r e , measured at the inter- 
mediate axis for the WFPC2 data and impose a typical 
model error of 3.5% (e.g., Paper V; the formal statistical 
errors are much smaller). 



We perform synthetic photometry using composite 
stellar population (CSP) models (see Section 5); mod- 
eling was performed separately for the HST photometry 
and the SDSS photometry, but we focus on the HST pho- 
tometry here. The several HST colors available for most 
systems allow us to constrain the restframe V-band lumi- 
nosity quite well, and the restframe B-band magnitude 
with slightly larger uncertainty. The CSP models are also 
evolved to z = to determine self-consistent evolved lu- 
minosities in the B and V bands; a more detailed look 
at the luminosity evolution will be presented in Paper 
X, but we note that the evolution inferred from the CSP 
models is in good agreement with other methods (Pa- 
per II, G09). Furthermore, we use the HST photometry 
and CSP models to perform synthetic photometry of the 
SDSS bands and compare these with the measured SDSS 
magnitudes for all systems with at least 3 HST bands (56 
systems in total) . Our comparison is with de-reddened de 
Vaucouleurs magnitudes (the deVMag magnitudes from 
the SDSS database corrected for Galactic extinction) in 
the bands g, r, i, z. We find that we are able to reproduce 
the SDSS photometry quite well, with typical random 
differences of only a few hundredths of a magnitude. To 
investigate the affect of the source light on SDSS photom- 
etry, we degrade our HST data to the resolution of the 
SDSS data (approximately 1"5) and perform the pho- 
tometric modeling on these data (this is done without 
adding additional noise; we therefore are only investi- 
gating the contribution of the background source in low- 
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resolution imaging and without explicit masks during the 
modeling). We find that we are able to recover the full- 
resolution magnitudes quite well; in general the sources 
are too faint to significantly affect the lens light photom- 
etry. 

We have not determined aperture colors from the HST 
photometry for the SLACS lenses (all of our HST pho- 
tometry is from full de Vaucouleurs model fits), but as 
a sanity check we use the SDSS photometric data to de- 
termine the restframe g — r aperture colors (here we use 
the S DSS modelMag magni tudes to determine the colors, 
e.g., iStoughton et al J 1 20 02). As is typical of early- type 
galaxies, we find a tight sequence with g — r ~ 0.8, and we 
plot a restframe color-magnitude diagram for the lenses 
in Figure [4j The absence of any significant tilt in the 
color-magnitude relation is due to the narrow effective ve- 
locity dispersion range of the SLACS lenses (comparable 
in width to the 'cons tant velocity dispersion' bins used 
by ot her authors, e.g.. lBernardi et aLll2005t [Graves et al.l 
2009); a slightly clearer trend is seen in a color-cr diagram 
(Figure EJ). 



0.95 




Figure 4. The restframe color-magnitude diagram for SLACS 
lenses. Black circles denote lenses at z < 0.15, blue crosses indicate 
lenses with redshift 0.15 < z < 0.2, green squares are lenses at 
0.2 < z < 0.25, and red diamonds are for lenses with z > 0.25. 
Typical early-type galaxies in the local universe have colors g — r ~ 
0.8. 



5. STELLAR MASSES AND STELLAR POPULATIONS 

One of the key benefits of having multi-color imaging 
is the ability to determine robust stellar mass estimates, 
particularly when near-infrared bands are included. The 
stellar masses of a subset of SLACS lenses have previ- 
ously been inv estigated using th e ground-based SDSS 
imaging (G09; IGrillo et alJ l20f)l . but the HST-based 
photometric data are substantially deeper, provide sig- 
nificantly better angular resolution, and extend into the 
near infrared for most objects. 

We have developed a novel code to perform a Bayesian 
exploration of the stellar populations of galaxies us- 
ing composite stellar po pulation models produce d with 
third-party codes (e.g.. iBruzual fc Charlotl [2003. here- 
after referred to as BC03). The code is described in de- 
tail in the Appendix, where we test the robustness of 
the code and the effects of parameter degeneracies and 
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Figure 5. The restframe color-tr diagram for SLACS lenses. The 
color /symbol coding is the same as in Figure [4] 

different filter combinations, but we briefly present the 
salient features here. The code takes a set of photometric 
data (magnitudes and errors through given bands) for a 
galaxy at a known redshift and creates a set of interpo- 
lation models for the synthetic magnitudes derived from 
a grid of CSP model spectra. Because the CSP models 
can be evaluated at arbitrary points within the bounds of 
the grid, the user can place arbitrary priors on the CSP 
model parameters. The parameter space is then explored 
using a Markov chain Monte Carlo (MCMC) routine, al- 
lowing a full determination of the posterior probability 
distribution function. 

In this paper, we restrict ourselves to using BC03 to 
create our CSP models, employing Chabrier or Salpeter 
IMFs and an exponential-tau model; the choice of models 
and IMF will be discussed elsewhere (Treu et al. 2009, 
in preparation). This star formation history model in- 
cludes five free parameters: the time that has passed 
since the beginning of the star formation (t), the char- 
acteristic time scale of the exponential burst (r) , metal- 
licity (Z), reddening (via the ry parameter of BC03), 
and the stellar mass (M*; see Table 5). We place a uni- 
form prior on t such that the star formation begins be- 
tween redshifts z = 1 and z = 5. An exponential prior 
is used for r with scale length I Gyr, and we use an 
extinction prior that is uniform in log[ry]. There is a 
known velocity di spersion-metal l icity c orrelation, and we 
use the sample of lGallazzi et al.l (|2005f ) to create a Gaus- 
sian prior on log[Z]; this is done by finding the mean and 
standard deviation o n log[Z] for all of the galaxies from 
IGallazzi et al.l (|2005f ) with similar redshift and velocity 
dispersion to the lensing galaxy. 

Table 5 

Stellar Mass Priors 



Parameter 


Type 


Range 


M» 


Log-Uniform 


10 9 to 10 13 M Q 


t 


Uniform 


t z= l to t z = 5 


T 


Exponential 


scale length of 1 Gyr 


T V 


Log-Uniform 


0.01 to 2.0 


z 


Log-Normal 


See text. 
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We apply this code to the HST photometry described 
in Section 4 and Table [3] The results of our analysis are 
listed in Table [H where we report the mean and stan- 
dard deviation of the marginalized posterior for the stel- 
lar mass for all systems with at least 2 HST bands; we 
also report the inferred stellar mass fraction within the 
Einstein radius (/*,Ein) f° r lenses with SIE models. Note 
that the stellar mass estimates do not include any priors 
from the lensing, and the inferred mass fractions within 
the Einstein radius may therefore exceed the total mass 
inferred from lensing (i.e., /*,Ein may be greater than 
1). We additionally provide the ratio of the fraction of 
light within the Einstein radius and the effective radius 
to allow easy conversion to the stellar fraction within the 
effective radius; note that this conversion implies an ex- 
trapolation of the lensing mass to larger radii under the 
assumption of an isothermal mass density profile. 

We have compared the stellar mass estimates derived 
using the HST photometry with stellar mass estimates 
from SDSS photometry. We first compare with the 
masses determined by the MPA/JHU group 10 . This sam- 
ple includes 68 of the systems with more than one band of 
HST imaging, and we find that the stellar mass estimates 
from the SDSS data are in very good agreement with our 
mass estimates using HST photometry (Figure H]) . We 
also use our stellar mass estimation code to determine 
masses from the SDSS photometry and find that we are 
able to reproduce the HST masses and the MPA/JHU 
masses. 
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Figure 6. Comparison between stellar masses determined us- 
ing the SDSS photometry and using the HST photometry. The 
red crosses are a comparison with the masses determined by the 
MPA/JHU group using SDSS photometry while the blue circles 
are a comparison with the stellar masses determined by G09 from 
SDSS photometry. Note that G09 tend to under-estimate the stel- 
lar masses of lower mass lenses compared to the MPA / JHU group 
and the masses from HST photometry. The red dotted line is a 
linear fit to the relation between MPA/JHU and HST masses, the 
blue dotted line is a linear fit to the relation between G09 and HST 
masses, and the black dotted line indicates the identity. 

10 Available at http://www.mpa-garching.mpg.de/SDSS/ 



G09 estimate the stellar masses of a subsample of 57 
of the SLACS lenses using four different CSP models, 
including a BC03 delayed exponential-tau model with 
metallicity set to solar metallicity and without modeling 
the effects of dust. Figure [6] shows the comparison be- 
tween the G09 stellar mass estimates and our stellar mass 
estimates for 55 of the systems investigated by G09 (we 
only have single band HST imaging for two of the systems 
in G09). We find a slight mass-dependent bias between 
the G09 stellar masses and those derived from our HST 
imaging; the G09 estimates for lower mass lenses are sys- 
tematically below the masses estimated from HST imag- 
ing. We note that galaxies with velocity dispersions simi- 
lar to the SLA CS lenses generally have super-solar metal- 
licities (e.g.. iGallazzi et al]|2005f ) and the G09 choice of 
solar metallicity may therefore be inappropriate; we find 
a range of metallicities, including some sub-solar, in our 
modeling. Differences in priors on the age-determining 
parameters (t and r) could also account for the offset 
seen between the G09 estimates of M* and the HST and 
MPA/JHU masses. 

6. THE PROPERTIES OF SLACS LENSES 

The final set of SLACS lenses spans a redshift range 
from z = 0.063 to z = 0.513, with a median redshift 
of Zmcd = 0.19. There is a strong correlation between 
the redshift and the mass of the lensing galaxy (as pa- 
rameterized by the Einstein radius rEi n ; see Figure [7]). 
This is largely the result of how SLACS lenses are se- 
lected; the SDSS spectra constitute a magnitude- limited 
survey 11 and the highest redshift objects are therefore 
more luminous (Figure [8|) and hence more massive. 




• • * « • 
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Figure 7. The redshift and Einstein radius, rEin, distributions for 
SLACS lenses. The median redshift is 0.19 and the median rEin 
is 4 kpc. There is a strong correlation between mass and redshift 
that results from the magnitude limit of the SDSS spectroscopic 
survey (see Figure [§}. 

The 85 SLACS lenses span approximately one decade 
in stellar and total mass, and are approximately evenly 
distributed across this range. The mean stellar mass-to- 
light ratio (M/L) in the V band is 3.2 (5.7) for a Chabrier 
(Salpeter) IMF; the average total M/L within half of r e 

11 Note, however, that the SDSS luminous red galaxy survey has 
a fainter magnitude limit as well as color criteria. 
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Figure 8. The redshift and magnitude distributions for SLACS 
lenses. The strong correlation between luminosity and redshift is 
the result of the selection function of the SLACS lenses, which arc 
drawn from the magnitude-limited SDSS spectroscopic survey. 

i r e/2) derived from lensing is 8.7 for either IMF, indi- 
cating that some dark matter must be present in both 
scenarios. We show the distribution of M* and the stellar 



mass fraction within r e / 2 , f, 



e/2' 



for the SLACS lenses 
in Figure [H We use r e / 2 because it is well matched to 
the Einstein radii of the lenses and therefore does not de- 
pend strongly on the profile of the lens model (e.g., Paper 
VII) . The lack of a strong correlation between these pa- 
rameters, and correlations between these parameters and 
the total mass and stellar velocity dispersion, will be ex- 
plored in more detail in Paper X. The mean stellar mass 
fraction within the Einstein radius of the lenses is 0.4 
(0.7) with a scatter of 0.1 (0.2) for a Chabrier (Salpeter) 
IMF. In principle, the combination of lensing and dy- 
namics can independently probe this fraction and will 
therefore discriminate between different IMF models; we 
leave this analysis for a forthcoming paper (Treu et al. 
2009, in preparation). 




10.5 11.0 11.5 12.0 

log[M./M s ] 

Figure 9. The distributions of the stellar mass, M», and stellar 
mass fraction within r e /2, f*,r e / 2 - There is no strong correlation 
between stellar mass fraction and stellar mass; the relationship 
between the stellar mass and other observables is explored in more 
detail in Paper X. 

As with previous SLACS papers (Paper II, Paper V, 



Paper VIII; also see lAugen I2008D . we want to test any 
biases of the SLACS lenses compared to the parent pop- 
ulation of SDSS galaxies from which they are drawn. In 
the context of this paper, this constitutes a comparison of 
the stellar masses of SLACS lenses to the stellar masses 
of SDSS galaxies with similar properties. We use the 
MPA/JHU SDSS stellar mass catalog to find stellar mass 
estimates for 100 SDSS galaxies which have a redshift 
within Sz/z < 0.1 of the lensing galaxy and a velocity dis- 
persion that is within 10% of the lens velocity dispersion 
(twelve systems did not have 100 comparison galaxies or 
did not have masses in the MPA/JHU catalog and are 
therefore excluded from this analysis) . The distributions 
of stellar masses for the SLACS lenses and for the com- 
parison sample is shown in Figure [TOl the SLACS lenses 
have stellar masses which are indistinguishable from their 
parent population. 
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Figure 10. The distributions of stellar masses for the SDSS 
comparison sample (black, open) and the SLACS lenses (green, 
hatched). The distributions are indistinguishable (the K-S p- value 
is ~ 0.95), strengthening our conclusion that SLACS lenses do not 
appear to be biased with respect to massive early-type galaxies 
with similar stellar velocity dispersions. 



It has been suggested that strong gravitational 
lenses are likely to be more concentrated in mass 
than non-lensing galaxie s (e.g.. Ivan de Ven et alj [20081 : 
Mandelbaum et al. 2008). This increased concentration 
may result from triaxiality, where the strong lens is pref- 
erentially aligned along the long axis. If this is the 
case for the SLACS lenses, and if the light follows the 
same distribution as the total mass, we would expect the 
SLACS galaxies to have smaller ellipticities than a com- 
parison sample of galaxies (preferential alignment along 
the long axis will, on average, decrease the ellipticity) . 
Note that the implications of triaxiali ty also depend on 
the number of lensed im ages formed (|Rozo et alj|2007t 
IMandelbaum et alj I2008D and two -image lenses, four- 
image lenses, and non-lenses would each yield different 
distributions of ellipticities if there is a significant bias 
due to triaxiality. However, Paper VIII has shown that 
the SLACS lenses are dominated by two-image lenses 
(approximately 85% are classified as doubles; SLACS 
does not appear to have a bias towards finding quad 
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lenses), and we therefore investigate the ellipticities of 
the full sample of lenses in comparison with non-lenses. 

We show the distribution of axis ratios q for SLACS 
lenses (for consistency we use the SDSS axis ratios, but 
we find that these are very close to the axis ratios inferred 
from the HST photometry) compared to the distribution 
of axis ratios for the comparison sample of galaxies in 
Figure [TTJ We do not find that the SLACS lenses are 
anomalously circular or elliptical; there is no evidence 
from the light distribution that SLACS lenses are more 
concentrated than the comparison sample, although the 
effect should be quite small and may b e unobservable 
with so few lenses (e.g. jRozo et al.ll2007l) . A decomposi- 
tion of the luminous and dark matter components would 
provide dark matter concentrations that could be com- 
pared with simulations, but that test is beyond the scope 
of this paper. 



4.0 




Axis ratio q 

Figure 11. The distributions of axis ratios for the SDSS compar- 
ison sample (black, open histogram) and the SLACS lenses (green, 
hashed histogram). The SLACS lenses do not appear to be more 
circular than the comparison sample (i.e., the green histogram is 
not skewed toward 1), indicating a lack of evidence for a preferred 
axis for the lenses. The two-sample K-S probability is ~ 0.13 for 
the two distributions, indicating that the two samples cannot be 
distinguished at greater than ~ 1.5<r confidence. 

7. DISCUSSION AND SUMMARY 

We have presented multi-band HST observations of the 
complete sample of SLACS gravitational lenses. These 
data have been used to upgrade 6 lenses from Paper V 
from possible to genuine (grade 'A') lens systems as well 
as confirming 9 newly investigated lens systems. We also 
present a novel fully-Bayesian stellar mass estimation 
code and have used this code with our HST photome- 
try to derive stellar masses for each of the grade 'A' lens 
systems. The new results of this work include: 

• The SLACS program has confirmed 15 new grade 
'A' lenses. SLACS has now discovered 85 confirmed 
lenses and has found 13 high-quality lens candi- 
dates, yielding a total of 98 likely lens systems. 
We have created SIE lens models with paramet- 
ric sources for 11 of the new lenses to allow for 
accurate determination of the mass within the Ein- 
stein radii of the lenses. Approximately 80% of the 



lenses have early type morphologies, while ~ 10% 
have spiral structure and the other ~ 10% have SO 
morphologies. 

• HST photometric colors have been used to deter- 
mine accurate and precise stellar masses for the 
SLACS lenses. Typical statistical errors on the 
mass estimates are 0.1 dex, and the masses de- 
termined using HST photometry are equivalent to 
those estimated from SDSS photometry. We note, 
however, that G09 find systematically lower masses 
for the less massive systems (~ 0.1 dex lower at 
10 11 M ). Assuming a Chabrier IMF, the mean 
stellar mass fraction within the Einstein radius of 
the lenses is 0.4 with a scatter of 0.1; the fraction 
increases to 0.7 with a scatter of 0.2 for a Salpeter 
IMF. The total mass-to-light ratio within half of 
the effective radius is 8.7. 

• There are no substantial differences between the 
stellar masses and/or ellipticities of SLACS lenses 
and a comparison sample of SDSS 'twin' galaxies. 
Thus, the SLACS sample appears to be, in effect, a 
velocity dispersion-selected sample of galaxies with 
the same properties as similarly-selected compari- 
son samples from SDSS: the SLACS lenses follow 
the same Fundamental Plane relationship (Paper 
II; Paper V II), reside i n the same environments 
(Paper VIIL lAugerl2008h . and have the same distri- 
butions of stellar mass and ellipticity as non-lensing 
galaxies (this work). 
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APPENDIX 

STELLAR POPULATION ANALYSIS CODE 

We present a novel Bayesian stellar population analysis code that allows the full posterior of the composite stellar 
population (CSP) model parameter space to be efficiently probed under the assumption of arbitrary priors, therefore 
providing error estimates for all quantities derived from these models (including evolved luminosities and so-called 
'k-corrections'). Additionally, we provide a routin e to perf o rm qu ick, accurate calculations of the Bayesian evidence 
using the nested sampling algorithm suggested by [Skilling (20Q3); this allows different model families (i.e. different 
choices of IMF) to be compared quantitatively. In this Appendix we describe in detail how the code works, we present 
the various tests we have used to ensure the stability of the code, and we discuss the effects of choosing different priors 
or different photometric bands when exploring stellar populations. 

Exploring the Stellar Population Parameter Space 

Conventional stellar population analyses tend to use a grid or discrete space of CSP models when investigating stellar 
populations. By cleverly choosing the axes of the grid or employing an a posteriori weighting scheme, informative 
(but discrete) priors can be placed on the various CSP parameters; however, it is also possible to inadvertently choose 
pathological priors. Our code circumvents the discrete priors by employing an interpolation scheme over a grid of 
CSP models that allows the models to be evaluated at arbitrary points within the bounds of the grid. For example, 
the analysis of SLACS lenses presented in Section 5 imposed a Gaussian prior on the metallicity of the CSP models 
informed by the velocity dispersion-metallicity correlation. The modeling is performed in a fully Bayesian framework 
and the user is therefore required to explicitly define priors for each parameter. 

The code operates in two stages. First, the user supplies a set of filters and a redshift for the object (future versions 
will likely allow the redshift to be a free parameter). The code then evaluates a model magnitude in each filter on a 
grid of (appropriately redshifted) spectra produced by CSP codes. These grids can have arbitrary dimensionality (e.g., 
any number of parameters among age, metallicity, etc.), although for practical computational purposes grids should 
generally be restricted to dimensionality < 6. A spline interpolation model of this grid of AB magnitudes is generated 
for each filter, and the set of models are then written to disk. The CSP models are undefined outside of this grid, 
and therefore care must be taken to make the grid sufficiently wide to allow for any physically interesting values of 
parameters. Furthermore, although the spline interpolation scheme and the smooth nature of most parameters allows 
for relatively sparsely sampled grids, the density of points should also be chosen to minimize interpolation errors. 

These interpolation models can be evaluated very quickly, allowing data to be efficiently analyzed via Markov chain 
Monte Carlo methods. The second stage of the code, then, is to define priors for each of the parameters in the model 
and the stellar mass and to explore the posterior of the parameter space given a catalog of photometric measurements. 
Note that while many codes use colors (i.e., the difference in magnitude between two filters) as the input data, our code 
uses the filter magnitudes and it is therefore possible to infer masses from single band observations; note, however, 
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that inferences from single band data are not particularly robust. Future versions of the code will likely add the ability 
to use fluxes to allow for the input of non-detections. 

The exploration of parameter space is undertaken in the usual Bayesian manner. The code uses a likelihood function 
that assumes Gaussian errors on the photometric magnitudes, 

filters , . . 

ln£(0) = - - 2^ — 2 + M 27rcr mJ; 

where 9 is the vector of model parameters, m; is the magnitude in filter i, modeli{6) is the CSP model magnitude 
evaluated at 0, and c^. is the observational error on the magnitude. The evaluation of the posterior, 

priors 

v$<xc(fi n p(0o 

i 

where P(0) are the priors, is handled by the Python package PyMC. The amount of time required to probe the 
parameter space depends on the number of data points and the number of components in the CSP model; a converged 
chain for a typical model with photometric data in four filters is generally obtained in less than one minute using a 
single processor on a modern workstation. 

Comparison with Other Methods 

Our code was origi nally tested agains t stellar masses determined by more conventional methods for galaxies in the 
GOODS -North field (iBuhdv et al.l l2005h. We accurately recover the stellar mass estimates (Figure 12) when the same 
priors as iBundv et al.l (|2005f ) are used ( i.e., implicit discrete priors on the metallicity and reddening) or when uniform 
priors bounded by the extrema of the IBundv et al.l (2005) priors are used. For these data and priors, the discrete 
and uniform priors lead to the same most likely value for the stellar mass; however, the uniform priors provide more 
accurate estimates of parameter uncertainties and covariances. We have also tested our code against the sample of 
SLACS 'twin' comparison galaxies discussed in Section 6 and find that we recover the MPA/JHU mass estimates to 
within 0.02 dex. 

Effects of Photometric Bands 

The stellar mass of galaxies is most robustly determined using near-infrared bands (e.g.. lBell &; de Jongll2001l but 
also see Maraston 2005; Conroy et al. 2009 for a discussion of potential systematics in the near infrared), and we now 
investigate the effects of increasing the number of filters for a given object. In Figure 13 we show the effects of including 
the H-band photometry on error estimates for all of the SLACS systems imaged by NICMOS. The red histogram 
shows the improvement gained when adding H-band data to BVI data, while the blue histogram demonstrates the 
improvement in precision when adding H-band data to a single optical color (i.e., two optical bands, either B-I or V-I). 
The inclusion of H-band data tends to decrease the errors by 20-30%, an improvement that partially results from the 
posterior of the stellar mass becoming more Gaussian when H-band data are used. 

Figures 14 and 15 show the effects of increasing the number of filters on the various parameter degeneracies for one 
object imaged in four bands, SDSSJ0037-0942. The metallicity constraint with just V and I photometry is weak and 
tends to reproduce the prior, whereas the addition of H-band data (Figure 14) narrows the inference and slightly shifts 
the distribution from the peak of the prior. The H-band data also eliminate models with significant dust, although 
this is not peculiar to the H filter since going from two bands to three bands tends to improve the inference on 
dust reddening in general. The exclusion of dusty models shifts the peak of the posterior of the stellar mass for this 
object, decreasing the inferred mass by approximately 0.1 dex. The addition of the B-filter data further constrains the 
reddening by dust, and the B band also eliminates models with extended star formation histories, again improving the 
precision of the stellar mass estimate (Figure 15). 

Effects of Priors 

We have seen in Section 5 that the choice of priors can noticeably affect the inferred stellar mass estimates. We have 
elected to use informative priors for all of our model parameters except the stellar mass, for which we use a broad 
uniform prior (note that the uniform prior on the log of Ty adds more weight to low optical depths, therefore encoding 
our belief that early- type galaxies are unlikely to be significantly affected by dust reddening). One could use even 
more constraining priors (by using the correlation between velocity dispersion and t suggested by I Thomas "et~aT1[200a 
for example); however, we now investigate the affects of using less informative priors. 

We begin by investigating the prior on the metallicity derived from other SDSS early-type galaxies (Gallazzi ct al. 
2005). We replace this with a uniform prior on the log of the metallicity, ranging from Z = 0.0001 to Z = 0.05. Figure 
16 shows the effects of the prior on the inferred metallicity; the prior affects the inference for all systems, but most 
notably for the systems imaged in only two bands (shown as blue diamonds in Figure 16). The corresponding changes 
in the stellar mass are less dramatic (Figure 17); M* is robustly determined when three or four photometric bands are 
used. However, the strong metallicity prior is particularly important for the systems with only two bands of data, as 
these systems tend to underestimate the stellar mass by ~ 0.1 dex without the prior. 
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Figure 12. Dis tribution of o f fsets i n stellar mass estimates for GOODS galaxies as determined by our new stellar mass estimation code and 
the estimates of Bundy ct al. (2005), assuming a Chabrier IMF. There is a negligible bias, AM= —0.03, and the width of the distribution 
is comparable to the errors on individual mass estimates, a 0.15. 



We also investigate the affects of using 'ignorant' priors by loosening the priors on t, the time since the beginning 
of the star formation, and t, the characteristic time scale of the exponential star formation rate model. We replace 
these with uniform priors, from 0.6 Gyr to the age of the Universe for t and from 0.04 Gyr to 5.1 Gyr for r; we also 
use the uninformed prior on the metallicity described above. The comparison between informative and ignorant priors 
in shown if Figure 18. The effect of the metallicity prior is evident, but we also find that removing the priors on the 
age of the stellar population causes the systems observed in only three bands to also shift. The uninformative priors 
allow for younger stellar populations to be present, thereby underestimating the stellar mass compared to the priors 
that require somewhat older stellar populations. 

The clear effects of the priors underscore the i mportance of choosin g reasonable priors when inferring stellar masses. 
Our set of ignorant priors mimics the priors in iBundv et al.1 <|2005f ). but we note that those authors investigated a 
heterogeneous set of galaxies over a range of redshifts, with four bands extending redward only to the observed-frame 
z-band, and the broad priors are therefore appropriate (at least to the extent that using a single set of priors is 
appropriate). The informed priors used in our analysis of the SLACS lenses are well- motivated, as the SLACS lenses 
are generally massive early-type galaxies at moderate redshifts (z ~ 0.2), which typically have quiescent, older stellar 
populations. 
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Figure 13. Demonstration of the improvement in errors gained when near infrared bands are used to estimate stellar masses. Actm, is the 
difference between <r pt, the error on the stellar masses of SLACS lenses when only including optical data, and <T a n, the error when the full 
photometry, including the F160W magnitude, is used. The red histogram is the improvement for lenses with imaging in four bands (three 
of which are optical) and the blue histogram shows the improvement for lenses with three bands of imaging (two of which are optical) ; the 
purple region is the overlap of the two distributions. 
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Figure 15. Joint marginalized distributions for the CSP parameters comparing B-V-I 
diagonal) to constraints from V-I photometry (lower-left corner, blue on the diagonals); 
photometry shown in Figure 14. 
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Figure 16. Comparison between the inferred metallicity of SLACS lenses assuming an informative or ignorant prior, where the informative 
prior is described in Section 5. The dotted line is the identity (note the different axis scales) while the blue diamonds are systems with two 
bands of imaging, the red crosses were observed in three bands, and the black circles have four bands. Solar metallicitv is -1.7. 
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Figure 17. Comparison between the inferred stellar masses of SLACS lenses assuming an informative or ignorant prior on the metallicity. 
The dotted line is the identity while the blue diamonds are systems with two bands of imaging, the red crosses were observed in three 
bands, and the black circles have four bands. 
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Figure 18. Comparison between the inferred stellar masses of SLACS lenses assuming strong or weak priors on the metallicity and star 
formation history. The dotted line is the identity while the blue diamonds are systems with two bands of imaging, the red crosses were 
observed in three bands, and the black circles have four bands. 
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Note. — Mb and My arc rcst-framc Johnson B and V AB magnitudes inferred from our CSP modeling. Mb,o and My.o arc the same, except the CSP models have been evolved to z 
(the ages have been increased by the amount of cosmic time between the redshift of the galaxy and z — 0, and the CSP models are re-evaulated with these new ages). 



